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SUMMARY 


Calibration data for the two-dimensional test section of the Langley 
0.3-Meter Transonic Cryogenic Tunnel are used to develop a Mach number-Reynolds 
number correlation for the fan pressure ratio in terms of test-section condi- 
tions. The relationship developed provides an excellent correlation of these 
data with the tunnel test section empty over a Mach number range from 0.30 to 
0.85 and a range of Reynolds number values which vary by a factor of 37. Well- 
known engineering approximations are used to derive an equation which is func- 
tionally analogous to the correlation deduced from these data. The success of 
the correlation is partially dependent on correctly calculating the fan pres- 
sure ratio from the static pressure measurements made on either side of the fan. 
The procedure for this calculation is described in the appendix. 

Although initially developed to simplify automatic tunnel control schemes, 
the correlation also provides (for this particular configuration) a successful 
Reynolds-number scaling law valid over the factor of 37 in Reynolds number. 

Also, a loss coefficient is formed which is independent of Mach number and 
Reynolds number and is a function of tunnel geometry only. This geometric loss 
coefficient should be useful as a measure of tunnel efficiency for operationally 
similar tunnels. 


INTRODUCTION 

The development of the cryogenic wind-tunnel concept (refs. 1 and 2) has 
greatly expanded the subsonic and transonic test capabilities available to the 
researcher. In addition to providing a large increase in Reynolds number, the 
ability to independently vary pressure, Mach number, and temperature makes it 
possible to perform the highly desirable research task of separating the aero- 
elastic, compressible, and viscous effects for the aerodynamic parameters being 
measured. For fan-driven tunnels at equivalent Reynolds number, those operated 
cryogenically require less power at the drive motor than conventional ambient 
tunnels. However, the power used in the production of the liquid nitrogen 
required for cooling tends to offset the savings in drive power. 

To utilize the cryogenic testing technique in an optimal manner, computer- 
ized tunnel controls are necessary to minimize the quantity, and thus the cost, 
of the liquid nitrogen used as coolant and also to reduce the time required per 
data point (ref. 3) . The development of tunnel control computer software 
requires knowledge of the functional relationship of the tunnel drive-fan power 
to Mach number, pressure, Reynolds number, and temperature. Any anomalies in 
the nature of the tunnel flow processes, particularly as a function of the wide 
variation in Reynolds number, must also be understood. This paper provides some 
insight into the control problem and the tunnel flow processes through a corre- 
lation of the fan pressure ratio for the Langley 0.3-Meter Transonic Cryogenic 
Tunnel with Mach number and Reynolds number. 



SYMBOLS 


A cross-sectional area open to flow 

k loss coefficient 

k Q loss coefficient independent of Reynolds number 

M Mach number 

• 

m mass flow 

n inverse of Reynolds number exponent 

p pressure, atm (1 atm = 101.3 kPa) 

R Reynolds number based on test-section hydraulic diameter 

R x local Reynolds number based on length 

^ gas constant 

r fan pressure ratio, Pt,2/Pt,l 

T temperature 

V velocity 

Z compressibility factor, p/p^R 

8 cooling capacity 

Y ratio of specific heats 

h efficiency factor 

p density 

Subscripts: 

1 upstream of fan 

2 downstream of fan 

f at fan 

s static conditions 

t total conditions 
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WIND TUNNEL 


The Langley 0.3-Meter Transonic Cryogenic Tunnel (TCT) is a single-return, 
fan-driven pressure tunnel with interchangeable test sections, including a 
20- by 60-cra two-dimensional test section which was in place during the present 
tests. Sketches of the 0.3-m TCT and two-dimensional test section are shown in 
figure 1 . The tunnel flow is maintained by a fixed-geometry single-stage fan 
driven by a 2.2-MW water-cooled synchronous electric motor. The motor, which 
is external to the tunnel, is powered from a variable-frequency source at fan 
speeds from 600 rpm to 5600 rpm. Further descriptions of the 0.3-m TCT may be 
found in references 3 and 4. 

A theoretical operating envelope for the 0.3-m TCT at a test-section Mach 
number of 0.85 is shown in figure 2 in terms of power delivered to the stream 
by the fan as a function of Reynolds number for various values of total pres- 
sure and total temperature. The equation for fan power in cryogenic nitrogen 
may be found in reference 5. To obtain the wide range of temperatures shown in 
figure 2, liquid nitrogen (LN2) is sprayed directly into the stream to offset 
the heat added to the stream as a result of the work done by the drive fan and 
the wall heat flux into the stream through the tunnel wall insulation. As a 
consequence of this direct method of cooling, gaseous nitrogen, rather than air, 
is the working fluid. An exhaust system which bleeds gaseous nitrogen from the 
low-speed end of the tunnel (fig. 1(b)) serves as the tunnel total pressure 
control. 


EXPERIMENTAL DATA 

The pressure ratio across the fan was computed from static pressure mea- 
surements taken from wall pressure taps upstream and downstream of the fan. 

The pressure measurements were taken during an initial calibration of the two- 
dimensional test section, which employs slotted top and bottom walls with 
5 percent open area. The test section was empty except for the calibration 
rake. Data were taken for selected combinations of three values of total pres- 
sure, three values of total temperature, and nine values of test-section Mach 
number. Total pressure was varied from 1.23 to 4.80 atm, total temperature from 
105 K to 301 K, and test-section Mach number from 0.30 to 0.87. The resulting 
test-section Reynolds number varied by a factor of 37 from 0.756 x 10® to 
28.0 x 10 6 per meter (2.48 x 10® to 91.8 x 10® per foot). 


DISCUSSION 
Data Analysis 

The fan pressure ratio may be expressed either as the static or the total 
pressure rise across the fan. For this analysis, the total pressures on each 
side of the fan were derived from the static pressure measurements as discussed 
in the appendix. The value of the total pressure ratio across the fan at each 
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calibration point is presented in figure 3 as a function of Reynolds number for 
the various nominal values of test-section Mach number. Inspection of these 
curves reveals, with one exception, a decreasing fan pressure ratio with 
increasing Reynolds number at a given Mach number as would be expected because 
of reduced viscous losses around the tunnel circuit. The exception noted occurs 
at M = 0.87. These data represent the tunnel operating in a choked condition, 
i.e., where an increase in fan speed and drive power produces no increase in 
test-section Mach number. In this condition, the pressure ratio is not related 
to Reynolds number, but rather to the fan-speed setting at which the tunnel 
operator judged the tunnel to be choked. Therefore, the fairing for the 
M = 0.87 data is dashed, since the individual data points cannot be considered 
to follow any logical trend. Also, solid symbols are used in figure 3 as well 
as several other figures to identify those points obtained under choked condi- 
tions. The data obtained under choked conditions are not considered in the 
ensuing calculations and discussion. 

The pressure ratio across the fan represents the sum of the tunnel circuit 
pressure losses. The results in figure 3 illustrate the normally expected trend 
of increases in tunnel efficiency (decreasing fan pressure ratio) with increas- 
ing Reynolds number. Data fairings in this figure also indicate that the reduc- 
tion in fan pressure ratio with increasing Reynolds number is gradual and well 
behaved; this result would be expected from decreasing viscous losses rather 
than more rapid changes such as might occur, for example, with improved diffuser 
performance due to a sudden reduction in the amount of diffuser separation. In 
order to more fully understand and predict the operating characteristics of the 
tunnel and to develop automatic tunnel control schemes, it is desirable to 
express the fan pressure ratio information in terms of the tunnel variables in 
as concise a form as possible. 

Polynomial curve fits of the data shown in figure 3 are one means of 
relating the fan pressure ratio to the other tunnel variables. However, because 
of the broad range of operating conditions, an awkward number of coefficients 
are required to express the data in equation form to the desired degree of accu- 
racy. Also, a simple polynomial curve-fit scheme which provides an acceptable 
match to the data over the entire range of variables was not found. As a result 
of these problems, an improved correlation scheme was sought. 

The data in figure 3 indicate that the fan pressure ratio is a function 
of Mach number and Reynolds number. The relative effects of these two vari- 
ables are shown in figure 4, where the data are plotted as a function of Mach 
number squared. As this figure indicates, the fan pressure ratio is approxi- 
mately linear with Mach number squared, and the spread in the data is orderly 
with Reynolds number. A simple power-law function for Reynolds number was 
assumed, and the exponent was varied until an optimum collapse of the data 
spread was obtained. As illustrated by the data presented in figure 5 (except 
for the data obtained under choked conditions) , this Mach number -Reynolds num- 
ber function gives an excellent correlation of the experimental data. For con- 
venience and for referencing of the results to test-section conditions, the 
Reynolds number was based on the test-section hydraulic diameter (test-section 
area multiplied by 4 and divided by test-section perimeter) . 
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Theoretical Considerations 

Although the relationship of fan pressure ratio to Mach number and Reynolds 
number described in the preceding paragraph was deduced from inspection of the 
data, a similar equation may be derived from conventional engineering relation- 
ships. As noted in reference 6, the pressure loss coefficient k may be writ- 

Ap A 

ten as k = , which is analogous to the drag coefficient in aerodynamic 

- pV 2 A 
2 

calculations, if A is the flow area in the region of the fan, and both Ap 
1 - 

and - pv z are evaluated at this region, then k may be rewritten as 


2(P2 “ Pi) 2 /P2 \ 2 

k = = 1 I = (r s - 1) 

YPi 2 YMjAPi / YMt 2 


where r s is the static pressure rise across the fan. This equation may now 

y 

be written in terms of r s ; that is, r s = k - M-| 2 + 1. 


Another engineering relationship, described in reference 7, expresses the 
friction loss coefficient as a product of an inverse function of a fractional 
power of the length Reynolds number and a constant of proportionality. An 
equivalent form of this expression for the present case might be k = koR -1 / n 
where k Q is the constant of proportionality. Reference 7 also presents values 
of the exponent n as a function of length Reynolds number and shows that the 
value of n increases with increasing Reynolds number. This can be interpreted 
to mean that as the Reynolds number increases, the amount of reduction in the 
loss coefficient decreases. Figure 15 of reference 7 shows that for an increase 
in length Reynolds number for 10^ to 10®, the value of n increases from 3 
to 9. In the present experiment, assignment of a length Reynolds number to the 
entire tunnel circuit is difficult. However, at a nominal test-section Mach 
number of 0.80, a total pressure of 3 atm, a total temperature of 200 K, and 
assuming a virtual origin of the flow at the tunnel screens, the length Reynolds 
number on the test-section wall is approximately 0.25 x 10®. This value basi- 
cally represents the length of the contraction section plus part of the test 
section and accounts for approximately 1/6 of the tunnel circuit length and a 
corresponding part of the circuit losses. It may be equally valid to assume 
the boundary-layer origin and termination at the fan, which increases the length 
by a factor of roughly 6. Thus, the length Reynolds numbers to be considered 
in choosing a value for n are on the order of 10®. The experimentally derived 
value for -1/n (-0.096; fig. 5) for the 0.3-m TCT circuit losses falls between 

-1/10 and -1/11, whereas reference 7 would predict a value of -1/8 to -1/9. 


5 



However, reference 7 presents results for the skin friction drag of flat-plate 
boundary layers and is probably not appropriate for estimating the effects of 
the more complicated tunnel circuit geometry. Reference 8 presents drag data 
for airships at length Reynolds number on the order of 1 0 9 which vary approxi- 
mately as R x -1 / 10 ; these drag data are also quoted to vary with values of n 
as high as Rx -1 /^ (fig. l and page 2 of ch. 14 in ref. 8). 


In light of the foregoing discussion, representation of the drag by a 
Reynolds number function of the type used here is not without precedent, and 
the value of the Reynolds number exponent is typical of other published data 
in this length Reynolds number regime. Thus the expression for fan pressure 


ratio might have been expected to have the form r 



k = k 0 R -1 /n f and the fan pressure ratio can be written r = 


+ 1 where 

^ koR-Vn m 2 


+ 1 


where n «» 10. Note that k must be a function of geometry and Reynolds num- 
ber, but k Q is a function of geometry only and thus represents a loss coeffi- 
cient which is independent of Reynolds number . 


The best fit to the experimental data as shown in figure 5 is 
r = 0.8205 m 2 R -0 *® 9 ® + 1.001 which is functionally analogous to the expres- 
sion derived in the preceding paragraph. The second term can also be written 
1 + 0.001, implying a small pressure-loss term independent of the Mach number- 
Reynolds number function, or a small systematic error in the data system. 

(Note that r may also be expressed as Ap t f/p t + 1 . ) This small pressure- 
loss term is difficult to resolve because the bookkeeping on pressure losses 
around the circuit is complicated by the cooling method, which in the case 
of the 0.3-m TCT, consists of spraying liquid nitrogen directly into the tun- 
nel circuit. An initial loss results from accelerating the injected liquid 
droplets, and a total pressure gain results from droplet evaporation. The 
aerothermodynamics of a similar process are discussed in reference 9. Prelim- 
inary estimates based on the information provided by reference 9 indicate a 
net gain in total pressure. However, there is also a pressure loss involved 
in exhausting gaseous nitrogen later in the tunnel circuit (the nitrogen 
exhausted is equal to the mass flow of the injected liquid nitrogen) . The small 
size of the independent term (0.001) indicates, at least to the first order, 
that the gains and losses in total pressure due to the liquid-nitrogen cooling 
and exhaust system are largely self- canceling. Additional measurements in the 
0.3-m TCT would be necessary to further resolve these second-order effects. 


Application 

A cryogenic tunnel of the type described in this paper is well suited to 
determining the value of tunnel-circuit loss coefficient, since the tunnel has 
a large Reynolds number range. The fan pressure ratio equation from figure 5 
mhy be used to serve several purposes. It allows liquid-nitrogen requirements 
and power usage to be linked in a relatively simple manner to conditions in 
the test section. Figure 6 depicts an example of a tunnel performance map pre- 
pared using the fan pressure rise correlation; the power delivered by the fan 
to overcome the tunnel pressure losses is plotted as a function of the fan 
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pressure ratio for a range of test conditions. This simplifies the scheduling 
of test points for minimum overall cost. Also, incorporation of the corre- 
lation in computerized tunnel control schemes improves efficiency and should 
lower testing costs. This is particularly true for anticipatory control algo- 
rithms which seek to arrive at a set point in the most efficient manner. 

In a broader sense, the coefficient of the Mach number-Reynolds number 
term, which is also the slope of the straight-line curve fit in figure 5, 
serves as a measure of the tunnel circuit efficiency; the larger the coeffi- 
cient, the less efficient the tunnel. Thus, the coefficient would be expected 
to vary with any change in tunnel losses such as drag due to models in the test 
section. These variations are expected to be second order in nature and, with 
increasing tunnel operation experience, should be predictable in terms of model 
type and attitude. Note that the coefficient is a function of geometry only 
and is independent of Reynolds number or Mach number. Thus, an expression of 
this type should serve to determine relative efficiency of comparable types of 
wind tunnels. 

An additional observation is that, in the absence of shock losses, the fan 
pressure ratio correlation serves as a Reynolds number scaling law for this 
particular wind-tunnel configuration. 


CONCLUDING REMARKS 

Calibration data for the two-dimensional test section for the Langley 
0.3-Meter Transonic Cryogenic Tunnel have been used to develop a Mach number- 
Reynolds number correlation for the fan pressure ratio in terms of test-section 
conditions. It has been shown that well-established engineering relationships 
can be combined to form an equation which is functionally analogous to the cor- 
relation. Additionally, a geometric loss coefficient which is independent of 
Reynolds number or Mach number can be determined. Present and anticipated uses 
of this concept include improvement of tunnel control schemes, comparison of 
efficiencies for operationally similar wind tunnels, prediction of tunnel test 
conditions and associated energy usage, and determination of Reynolds number 
scaling laws for similar fluid-flow systems. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
October 2, 1980 
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APPENDIX 


CALCULATION OF FAN TOTAL PRESSURE RATIO 
FROM MEASURED STATIC PRESSURES 

During tunnel calibration tests of the Langley 0.30-Meter Transonic Cryro- 
genic Tunnel, static pressures were measured at stations upstream and downstream 
of the fan in addition to the measurements which were made to determine the 
test-section flow parameters. The cross-sectional areas at the static-pressure 
measurement stations were determined from design drawings. This information is 
sufficient, in the region of the fan, for calculating the total pressure ratio 
across the fan if it is assumed that there is negligible heat transfer at the 
tunnel walls and that the stagnation temperature at the station downstream of 
the fan is the same as that at the test section. The stagnation temperature may 
not be constant because of the exhausting of gaseous nitrogen that takes place 
in between these two stations, but the total pressure ratio calculation is not 
sensitive to this assumption. 

To calculate the total pressure ratio, initial calculations are made with 
the assumption that liquid nitrogen is not being added in the circuit; thus, 
the mass flow at the fan is the same as the mass flow in the test section. The 
following equations are iteratively solved for the downstream total pressure, 
Mach number, and static temperature: 



T t ,2 = T 2 ^l + ^ M 2 2 ] 

rr 

P2 m 2 a 2 v ~ 

\l Z t 

The term z t is a real gas correction described in reference 5. The ideal 
diatomic gas value of 1.4 is used throughout this report for the ratio of 
specific heats. 

The upstream flow condition can be solved with the assumption of fan 
efficiency factor. 
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APPENDIX 


At' 



where At' is the ideal temperature rise across the fan assuming isentropic 
compression, and At is the actual temperature rise. The ideal temperature 
ratio across the fan is related to the pressure ratio by 

T t,l \Pt,!/ 


&i /Pt,2\ <Y ~ 1>A 
T t ,l \Pt,l/ 

where At^ - " T t,l 

Substituting for At^, 


At 


nf 


^(Y-D/T 

J>t,l/ 


1 


and 


T t,l 


T t, 2 

1 Fpt,2\ (Y ~ 1)/Y 

^fKPt,!/ 


+ 1 


This equation along with the equations 


APPENDIX 


Pt,l " 


Pi 

Y-l __ 2 \-Y/(Y-D 


1 + Mi 

2 



are iteratively solved for Pt,lf T t f l/ Mi, 911(3 T 1 • 

The stagnation conditions on both sides of the fan have been determined, 
but the addition of mass flow through the fan due to the injection of LN 2 
for cooling has been neglected. The fan power can now be calculated from 


Power 



*<Tt f 2 


T t ,l)Zt 


The mass flow of LN 2 
is calculated from 


Power 


m LN 2 



necessary to remove the energy being input by the fan 


where 3 is the cooling capacity of LN 2 as given in reference 10. This mass 
flow is added to the test-section mass flow, and the whole series of iterative 
calculations are repeated with this additional mass flow through the fan. This 
procedure is repeated until there is no appreciable change in the mass flow 
going through the fan, and in general, takes only two or three iterations. 

The calculations in this paper used a fan efficiency of 0.89. Calculations 
using other fan efficiencies indicate that the total-pressure ratio across the 
fan is insensitive to fan efficiency. 
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(a) Two-dimensional test section in place. 


Figure 1.- Langley 0.3-Meter Transonic Cryogenic Tunnel. 
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Figure 1.- Concluded. 
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Figure 2.- Fan power requirements as a function of Reynolds number for a test-section Mach 
number of 0.85. (Reynolds number based on test-section hydraulic diameter.) 
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Figure 6.- Fan power as a function of fan pressure ratio for typical 
ambient and cryogenic operating points, assuming a fan efficiency 
of 0.89. 
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